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Abstract. We consider a two-dimensional Fermi-Pasta-Ulam (FPU) lattice with 
hexagonal symmetry. Using asymptotic methods based on small amplitude ansatz, 
at third order we obtain a reduction to a cubic nonlinear Schrodinger equation (NLS) 
for the breather envelope. However, this does not support stable soliton solutions, so 
we pursue a higher-order analysis yielding a generalised NLS, which includes known 
stabilising terms. We present numerical results which suggest that long-lived stationary 
and moving breathers are supported by the lattice. We find breather solutions which 
move in an arbitrary direction, an ellipticity criterion for the wavenumbers of the 
^ ' carrier wave, asymptotic estimates for the breather energy, and a minimum threshold 

energy below which breathers cannot be found. This energy threshold is maximised for 
stationary breathers, and becomes vanishingly small near the boundary of the elliptic 
' domain where breathers attain a maximum speed. Several of the results obtained are 

similar to those obtained for the square FPU lattice (Butt & Wattis, J Phys A, 39, 
4955, (2006)), though we find that the square and hexagonal lattices exhibit different 
properties in regard to the generation of harmonics, and the isotropy of the generalised 
NLS equation. 
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1. Introduction 

Discrete breatliers are time-periodic and spatially localised exact solutions of 
translationally invariant nonlinear lattices. For a brief review of some general properties 
of breathers in higher-dimensional systems, see our earlier work [2I. In particular, it is 
known that while some fundamental properties such as the existence of breathers are 
not affected by lattice dimension (see Flach et al and Mackay and Aubry other 
properties are affected profoundly, for instance, the energy properties of breathers (see 



Flach et al [Sj, Kastner and Weinstein 19||). In this paper we investigate how the 
symmetry of the lattice influences the properties of discrete breathers found therein. 

In Hamiltonian systems, stationary breathers occur in one-parameter families. For 
a certain class of Hamiltonian systems, a critical spatial dimension dc exists such that 
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for systems with d > dc, there exists a positive lower bound on the energy of a breather 
family, and the breather energies do not approach zero even as the amplitude tends 
to zero. For lattices with dimension d < d^, there is no positive lower bound on the 
energy of breathers. In other words, the energy of a family of breathers goes to zero 
with amplitude, and breathers of arbitrarily small energy can be found. The critical 
dimension d^ is typically two. A small amplitude expansion yields the NLS reduction 
to iFt + AF + = which has a critical dimension of dc = 2/cr, with blow-up 

in NLS occurring when k > and d > dc- For typical lattice potentials, a = 1 again 
confirming the critical dimension of dc = 2. 

Marin, Eilbeck and Russell have performed extensive numerical investigations of 
breather dynamics in two-dimensional lattices [isl, 0, Their results suggest 

that moving breather modes exist (or are at least extremely long-lived), and that 
the lattice exhibits a strong directional preference whereby breathers can only move 
along symmetries of the lattice, and in no other direction. Such quasi-one-dimensional 
behaviour is also observed in higher-dimensional lattices. 

The work in this paper follows on from our earlier study (3| of breathers in a 
two-dimensional lattice with square rotational symmetry; hereupon referred to as the 
"square" lattice. This lattice has C4 symmetry, by which we mean that a rotation through 
any multiple of 27r/4 = 7r/2 about an axis perpendicular to the lattice plane through 
a lattice site maps the lattice onto itself. In using the semi-discrete multiple-scale 
method (see Remoissenet [l8|), we determined an approximate form for (as well as the 
properties of) small amplitude breathers in a two-dimensional square Fermi-Pasta-Ulam 
(FPU) lattice 0|. We found a third-order analysis to be inadequate, since the partial 
differential equation obtained at this order describing the breather envelope exhibits 
blow-up. To overcome this, we incorporated higher-order effects in the model, thereby 
obtaining a modified partial differential equation which includes known stabilising terms. 
From this, we determined regions of parameter space where breather solutions are 
expected. Numerical simulations supported the results of our analysis, and suggested 



that, in contrast to the two-component lattices studied by Marin et al [l5|, ll6|, ll7|| (that 
is, with two degrees of freedom at each lattice site), there is no restriction upon the 
permitted directions of travel within the one-component square FPU lattice. We also 
found asymptotic estimates for the breather energy which confirmed the existence of a 
minimum threshold energy, in agreement with the work of Flach et al 

In this paper, we consider a hexagonal electrical transmission lattice (HETL). This 
two-dimensional network possesses Cg (or hexagonal) rotational symmetry. That is, 
rotation through any multiple of the angle 27r/6 = 7r/3 about a lattice site maps the 
lattice onto itself. The HETL is shown in Figure [T], pictured from a point vertically 
above the plane of the lattice. We note from Figure [T] that geometrically, the HETL is 
an arrangement of tessellating triangles (not hexagons). Nevertheless, the arrangement 
in Figured] is referred to as "hexagonal" rather than "triangular" since these descriptions 
refer to its symmetry properties and not to the geometrical shapes which comprise 
the array (not all authors follow this convention). One might expect the analysis for 
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Figure 1: The 2D hexagonal electrical transmission lattice (HETL). 

the hexagonal lattice to be more involved, since it is geometrically more complicated 
in having more links emanating from each node. However, the hexagonal symmetry 
results in greater isotropy and hence simpler equations than those obtained for the 
square lattice. 

We derive the equations of motion and demonstrate a Hamiltonian formalism in 
Section 12. 2[ In Sections 12.61 and 12.71 we present two cases for which the hexagonal FPU 
lattice equations can be reduced to a two-dimensional nonlinear Schrodinger (NLS) 
equation with cubic nonlinearity. We consider lattices with a symmetric interaction 
potential, in which case reduction to a cubic NLS equation can be performed for moving 
breathers. We find an ellipticity criterion for the wavenumbers of the carrier wave in 
Section [2.6. 1[ A reduction to the cubic NLS equation can also be carried out for lattices 
with an asymmetric potential, provided we consider only stationary breathers. 

As expected we find a minimum energy below which breathers cannot exist in 
the hexagonal FPU lattice. We find that the energy threshold is dependent upon the 
wavevector of the carrier wave. It is maximised for stationary breathers, and becomes 
arbitrarily small near the boundary of the elliptic domain. 

The cubic NLS equation admits only unstable Townes solitons. Hence, in Section [3l 
we extend our asymptotic analysis to higher order and find an isotropic generalised NLS 
equation which incorporates known stabilising terms. Our analytic work is supplemented 
by numerical simulations presented in Section [H which suggest that long-lived stationary 
and moving breather modes are supported by the system. In Section [5l we discuss the 
results obtained in this paper. 
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2. A two-dimensional hexagonal Fermi-Pasta-Ulam lattice 
2.1. Preliminaries 

Before we derive the equations governing the HETL, we describe our scheme for 
indexing the lattice nodes and our choice of basis vectors. We introduce a rectangular 
lattice with basis vectors B = {i' , j'}, where i' = i = [1, 0]^ and j' = [0, h]'^ (j being 
[0, 1]^), illustrated in Figure El We use only half of the (m, n) indices, namely, those for 
which the sum m + is even. We choose an origin with coordinates (0, 0); the position 
of the site (m, n) is mi' + nj'. In order for the hexagonal lattice to be regular, we specify 
h = V3. 



(m-l,n + l) (m + l,n + l) 



(m — 2, n) 









n) 











(m + 2,n) 



L 



(m — l,n — 1) (m + l,n — 1) 

Figure 2: Labelling of nodes in the HETL with basis B = {i' , j'}. 

At every node of the HETL lies a nonlinear capacitor (not shown in Figure [T]) , 
and between every node and each of its six nearest neighbours is a linear inductor. 
An enlarged view of the area surrounding the capacitor at (m, n) is shown in Figure [3l 
where the capacitor is visible. The variable „ denotes the voltage across the capacitor 
(m, n) and Qm,n denotes the total charge stored on the capacitor at (m, n). Also, Im,n, 
Jm,n and Kjji,n are the currents through the inductors immediately to the right of site 
{m,n) in the directions ej = [2,0]-^", ej = [1, — -\/3]-^ and = [1, -\/3]^ respectively, as 
illustrated in Figure [3l 



2.2. Derivation of model equations 

To derive the equations relating current, charge and voltage in the lattice we apply 
Kirchoff's law 



m+l,n— 1 <'Tn,n 



' m,n 



dt ' 
dJm,n 



Vr 



m+l,n+l ^m,n 



Vrr, 



dt 

dKm 



dt 



where the inductance L is constant. Conservation of charge gives 



'm— 2,n ^m,n 



^m— l,n+l Jm,n ~\~ -^m— l,n— 1 -^m,?^ 



dQr 



dt 



(2.1) 
(2.2) 
(2.3) 

(2.4) 
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m+l,n+l 




K-fi+l,)i-l 

Figure 3: Enlarged view of the HETL at site {m,n). 



Differentiating (12.4p witfi respect to time, and tlien using f l2.ip -( l2T3D to find Im-2,n, 
i^n+i and Kyn_i^n-i, we liave 



dt2 



(Vm+2,n — '2Vm,n + Vm-2,n) + (Kn+ l,n— 1 '2^m,n ~l~ l,n+l) 

+ (Kn+1 )• (2.5) 

Equation (12.50 may be written in the abbreviated form LQm,n = {S] + 5j + S]^)Vm,n, 
where the centred second-difference operators are defined by 

^m,n ~l~ -^Tn—2,ni 
'^■^m,n ^m— l,n+l) 
2^m,n ~l~ ^m— l,n— !• 



4 


= ^mH 


-2,n ~ 




= ^mH 


-l,n-l 


f)^ A 


= ^mH 


-l,n+l 



(2.6) 
(2.7) 
(2.8) 



Here, Ajn,n is an arbitrary quantity referenced by two indices; 6], 6j and Sj^ are centred 
second-difference operators in the directions of ei, ej and ek respectively. Since the 
voltage Vm,n is known in terms of the charge, we reformulate (12. 5p in terms of the 
single quantity Qm,n- We invert the capacitor's nonlinear charge- voltage relationship 
Q = VC{V) (see equations (2.9)-(2.14) of Q for details), to obtain 

V{Q) = {Q + aQ^ + bQ^ + cQ^ + dQ')/Co, (2.9) 
where Co = C(0). Hence, the HETL equations (12.51) can be written as 
d^O 



dt2 



{Sj + 5j + 5l) [Qm,n + aQl^^ + bQi^^ + cQl^^ + dQ^J , (2.10) 



where m, n G Z and, by rescaling the time variable, we set LCq = 1 without loss of 
generality. Thus we have shown that the equation governing charge in the HETL (|2.10l) 
is a two-dimensional analogue of the Fermi-Pasta-Ulam equation 

d^Q, 



dt^ 



j+i) 



2W'{Q,) + W'iQj 



(2.11) 
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which is a Hamiltonian system that can be derived from both 

^ = EK + ^(<^i+i-<^i)' and H = J2l{P,^,-P,Y + W{Q,),{2.12) 
j j 
where Qj = — (j^j. 

The lattice equations (12.101) can be derived from the Hamiltonian 

H ^ ^ 2i.Pm+2,n Pm,n) ~l~ 2^Pm+l,n—l Pm,n) 
m,n 

Pm,nr + ^{Qm,n), (2.13) 

where T{Qm,n) satisfies T'{Qm,n) = V{Qm,n) given in ([231), hence 

T(Q) = + |aQ3 + IbQ^ + IcQ' + IdQ'. (2.14) 

We describe potentials which satisfy T(— Q) = T{Q) (that is, a = c = 0) as 'symmetric'. 
The variables Pm,n and Qm,n are conjugate momenta and displacement variables of the 
system; and eliminating Pm,n from the equations 

^ = -(5/ + Sl)Pm.,n, ^ = -T'(g^,.), (2.15) 

yields f l2.10p . In Section 12.41 we derive expressions for the energy of breathers given the 
small amplitude solutions which are obtained in the next section. 

2.3. Asymptotic analysis 

We apply the method of multiple-scales to determine an approximate analytic form 
for small amplitude breather solutions of (12.101) . with slowly varying envelope. We 
introduce new variables defined by 

X = em, Y = ehn, t = et and T = eH; (2.16) 

note the presence of the scaling factor h in the definition of Y. We seek solutions of 
( j2T0ll of the form 

QmAt) = ee'^F{X, Y, r, T) + e''G^{X, F, r, T) + e-'e'^G^iX, F, r, T) + 

e'e''^G2{X, F, r, T) + e^H^{X, F, r, T) + e^e'^Hi{X, F, r, T) + 
s^e^'^^H^iX, F, r, T) + eh^'^fH^^iX, F, r, T) + e^e'^hiX, F, r, T) + 
e^e^'^^hiX, F, r, T) + e^e^'^^hiX, F, r, T) + e^e^'^^hiX, F, r, T) + 
eV^ Ji(X, F, r, T) + ■ ■ ■ + c.c, (2.17) 

where the phase ip of the carrier wave is given by km + Ihn + ci;t (once again noting 
the extra factor /i), and k = [k, and uj(k) are its wavevector and temporal frequency 
respectively. We substitute the ansatz (I2.17P into the lattice equations (12.101) and equate 
coefficients of each harmonic frequency at each order of e. After much simplification, 
this yields the following system of equations: 



0(ee''^): 

u'^F = 4sin2(A;)F + 4sin2 ('^±^) F + 4 sin^ ftzlh.] p, (2.18) 
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T 



ujFr = 2 sin /c [2 COS /c + cos(lh)]Fx + 2/icos/i;sin(//i)Fy, 

u'^G2 = [sm\2k) + sm\k + Ih) + sm\k - lh)]{G2 + aF^) 

= [4 cos(2A;) + 2 cos A; cos{lh)]Fxx + 2/i^ cos k cos(//i)Fyy 
— 4h sin k sin{lh)FxY 



(2.19) 
(2.20) 



126 



. 2/,\ • 2fk + lh\ . o(k-lh\ 
sin2(A;) + sm^ I ^ J + sm^ I — — J 

sm [kj + sm — - — + sm — - — 



[F(Go + Go) + FG2] 



(2.21) 



3„3i^/)^ 



0(£V): 



sin^(3/c) + sin^ 



3k + 3lh 



'3k - 3lh" 



sm 



H3 + 2aFG2 + 6F3) , (2.22) 



Gott — 6 Gqxx + 2/i^Goyy + 



6 IFI 



XX 



yy. 



(2.23) 



Though each equation plays a similar role to its counterpart in the square lattice, 
equations (12.18p -( !2.23p are more complicated. Equation f l2.18p is the dispersion relation 
for the system (I2.10p . Since we are interested only in solutions for which -F 7^ 0, (|2.18lj 
yields 



4sin^(/i;) + 4sin^ 



'k + Ih' 



+ 4 sin 



'k-lh' 



(2.24) 



which does not simplify significantly. From equation (I2.19p . we determine the velocity 
of the travelling wave F, finding 

F{X,Y,T,T) = F{Z,W,T), (2.25) 

where Z = X — ur and W = Y — vr, and the horizontal and vertical velocity components 
{u and v) are found to be 



-2 sin(A;) [2 cos(A;) + cos(//i)] 



-2h cos{k) sm{lh) 



u = anQ V = . [Z.zo] 

Equation (12.261) . along with (12.241) enables the elimination of terms involving G\ from 
(I2.2ip which are not shown. We denote the angle at which the envelope F propagates 
through the lattice by \E', which is measured from the direction of the basis vector ei 
to the line of travel and hence is given by tan~^(t;/u), which in turn depends upon the 
wavevector k. For both the cases that we consider (namely, symmetric and asymmetric 
interaction potentials), we find constraints upon the wavenumbers k and / which affect 
the velocity components u and v. By taking k = 7r/3 and I = n/h we find u = v = 0, 
which corresponds to a static breather; and by choosing k, I values which circle this 
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point, we find breathers which can propagate in any direction (0 < \1/ < 27r); in other 
words, our analysis suggests that there is no restriction upon the direction of travel 
through the lattice. 

Our aim is to reduce (12.211) to a nonlinear Schrodinger (NLS) equation for F. Before 
this can be done, the quantities Go and G2 in f l2.2ip must be found in terms of F. As 
for the square lattice, it is straightforward to determine G2 from the algebraic equation 
fl2.20p . However, the partial differential equation (|2.23l) for Go can be solved for two 
special cases only; namely, symmetric potentials, in which case the reduction of (12.211) 
to an NLS equation can be completed even for moving breathers, and also asymmetric 
potentials, provided we confine our attention to stationary breathers. These two cases 
are considered in Sections 12.61 and 12.71 respectively, where we also use our breather 
formulae to generate estimates for the breather energy. 

2.4- Asymptotic estimates for breather energy 

The HETL is a lossless network, meaning that the total electrical energy is 
conserved. This quantity is related to the Hamiltonian (12.131) hj E = H/Cq and so 
is given by 

^ ~ 5Z ^m,n. = ^ I" 2^ {jm,n + Jm,n + ^^m,n) • (2.27) 

m.n m,n ^0 

The electrical energy e^.n associated with a unit of the lattice is (see Figure [3]) 

T{Qm.,n) ,li(j2 ,j2 ,^2 ^ 28) 
^m,n ^ ^ 2 V m,7i ^ '^m,n ' ^^m,n)- ^^.^o; 

To derive a leading-order estimate for the electrical energy, defined by 

-^0 = ^ e^|n ~ ^ 1^ ^ (-^m,n + '^m,n + ^m,,ri) i (2.29) 
m,n m,n ^ 

we use leading-order expressions for each of the terms in the terms in the summand of 
( 12:271) . The first term is Q^„/(2Go); from ((231), it follows that Kn,n ~ <5m,n/Go, and 
so to leading order, (12.270 agrees with the linear approximation to the energy in the 
capacitor being QV/2. To find leading-order expressions for the currents /m,„, Jm,n and 
Km,m we use equations (12.ip - (12.31) . which imply 

Qm+2,n Qm,n ~ ) (2. 30) 

Qm+l,n~l Qm.n ) (2.31) 

dKrnn 

Qm+l,n+l Qm,n ; (2.32) 

where LGo = 1. The currents are determined by substituting the expression for the 
breather Qm,n into (12.301) - (|2.321) and then integrating with respect to time. We obtain 
leading-order estimates for the energy of moving breathers in systems with symmetric 
potentials Section 12.6.21 and in Section 12.7.11 the energies of stationary breathers with 
asymmetric potentials. 
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2.5. The dispersion relation for the HETL 

In this section, we analyse the dispersion relation f l2.24p for the system ( 12.10p . A 
contour plot of u against k and / is shown in Figure [H Since w is periodic in both k and 
/, with period 27r along the fc-direction, and 2n/h in the /-direction, we consider only k 
and I such that {k, I) e = [0, 27r] x [0, 2Tr/h]. 

The function uj is minimised, and assumes the value zero, at the centre of the 
circular patterns in Figure [H The minima are located at (0,0), (27r, 0), (27r, 27r//i), 
(0,27r//i) and {n^n/h). Plus signs ('+') mark the points in (/c, /)-space at which u is 
maximised and takes the value u; = 3. The maxima lie at the centres of the equilateral 
triangles, the wavevectors corresponding to these points are denoted ki, . . . , kg, where 

ki = [tt/S, n/hf, = [2n/3, Of, k, = [An/S, Of, 

k4 = [57r/3, 71/ hf, kg = [47r/3, 27r/ hf, kg = [27r/3, 271 /hf. ^ ' 

The arrangement of the points corresponding to these wavevectors in {k, /)-space reflects 
the hexagonal symmetry properties of the lattice (|2.10l) . It may be verified using 
(12 .26^ that the velocity components u and v are both zero for each of the wavevectors 
{ki,...,k6}. 




1 2 3 4 5 6 

k 

Figure 4: Contour plot of a;(k); uj attains its maximum value of 3 at the points 
ki, . . . , kg, and is minimised at the points marked 'o', where u; = 0. 



2. 6. Lattices with a symmetric potential 

In this section, we consider lattices with a symmetric interaction potential, where 
T((5) is even and T'i^Q) has odd symmetry. This corresponds to a = c = in (I2.10p . 
Since there are no even harmonics for vibrations controlled by symmetric potentials, 
it follows that Go and G2 are both zero. In (12.211) . the term Frr is eliminated using 
Frr = u^Fzz + 2uvFzw + v'^Fww which is derived from ^(TM . This leads to (IMTD 
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being reduced to an NLS equation 

2iujFT + — A cos(2/c) — 2 cos k cos{lh) Fzz + v'^ — 2/i^ cos k cos{lh) Fww 

+ \luv + Ah sin ksm{lh)]Fzw + ^huo'^\F\^F = (2.34) 

where the velocities u and v are given by (12.261) . By applying an appropriate change of 
variables, we remove the mixed derivative term from (I2.34p . and reduce the equation to 
a standard form. To simplify the appearance of subsequent expressions, we denote 
the coefficients of Fzz^ Fww and Fzw by Di = — 4cos(2A;) — 2 cos A; cos (//i), 
-D2 = f ^ — 2/i^ cos k cos{lh) and D3 = 2uv + Ah sin k sm{lh) respectively. A suitable 
transformation is thus 

. hZ . h{2D,W-D,Z) 

^ = -JW ^ = / , (2.35) 



which implies f l2.34p becomes 

2icjFT + 3V2F + 36cj2|F|2F = 0, (2.36) 

where the differential operator V^F = + is isotropic in the (^, 77) variables. An 
approximation to the Townes soliton solution of (12. 36^ is given by equation ( 1A6I) in the 
appendix. Substituting the resulting expression for F into ( 12.17^ yields a leading-order 
expression for the breather 

Qm,n{t) = 2ea cos[km + Ihn + (cj + e'^X)t] sech{f3r) + 0{e^), (2.37) 

where a and /5 are determined as described by equation (1A5I) in the appendix, with 
D = 3/2lj, B = ?>hu/2. Further, r = \/^^ + rj^ is found in terms of the physical discrete 
variables m and n by inverting the transformations (12. 35^ and reverting back to the 
variables Z and W , using 

Ah'^e'^{D2{m-utf + Di{hn-vtf - D^{m-ut){hn-vt)) 



2 f-2 I 2 

(2.38) 

The terms Di, £)2 and Da are known from (12.341) . the velocities u and v are given by 
(12:261) and u; is given in (12:241) . 

Determining the domain of ellipticity. We confine our attention to elliptic NLS 
equations. We seek to determine the region V of (/c, /)-parameter space (that is, the two- 
torus = [(0, 27r)] X [0, 27r/h]) where the NLS equation (l234l) is elliptic. By definition, 
this equation is elliptic when < 4D1D2, where Di, D2 and are given in ( 12.34^ . 
Whilst the region V cannot be specified explicitly, it is simple to find numerically, and 
is illustrated in Figure O Defining the function e{k,l) = 4:Di{k,l)-D2{k,l) — D^^kJY, 
we are concerned with the region where e(/c, /) > 0. The subdomains have been labelled 
{Pi, ... , Ve} in Figure U where V = ViU ...UVq. 

Again, the hexagonal symmetry properties of the HETL are reflected clearly in the 
function e{k,l), which has six maxima (at which e{k^l) = 36), each lying at the centre 
of one of the closed curves in Figure [5l The maxima of e(/c, /) coincide with the six 
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Figure 5: The domain P = Di U . . . U in which the NLS equation (|2.34p is elhptic. 

maxima of Lj{k, I) shown in Figure namely, at the wavevectors {ki, . . . , ke} in T^; 
e(fc, /) is minimised (e(/c, /) = —48) at the six midpoints of the line segments which 
connect adjacent maxima. 




Figure 6: Contour plots of Lo{k, I) (left) and speed {Vu'^ + v'^, right) for wavevectors 
k G Pi. Contours are for 2.82 < w < 3 in steps of 0.02 and speeds from 
zero to 0.6 in steps of 0.1. The central point corresponds to u; = 3 and 
zero speed, as one considers wave vectors nearer to the edge of Pi, the 
frequency uj decreases and the speed increases. 
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Figure E] shows a contour plot for the frequency u in Vi. the ellipticity constraint 
only permits breathers with a relatively high frequency, that is, with lj > 2.82. This 
constraint in turn implies that not all breather envelope velocities are attainable, only 
breathers with speeds upto about 0.7 lattice sites per second are permitted; the plot 
on the right of Figure [6] shows that only breathers with speeds upto about 0.3 sites 
per second can move in any direction. There is then an intermediate range of speeds, 
between 0.3 and 0.7 sites per second where breathers can only move in certain directions, 
these correspond to the lattice directions. The larger the speed, the more restricted is 
the direction of motion. 



2.6.2. Breather energy. To calculate the leading-order energy Eq of moving breathers 
in lattices with a symmetric potential we use f l2.37p which we write as 

Qm,n{t) ~ 2£:acos<l> sech(/3r), (2.39) 

where $ = km + Ihn + VLt is the phase of the carrier wave, {k,l) eV^VL = uj + e^X is the 
breather frequency including the first correction term, = km + Ihn + cot is only the 
leading-order expression). A parameterises the breather amplitude, a = a(A), /3 = /3(A) 
and are as described in equations (12.371) and (I2.38p . 

We now find expressions for the currents Im,n, Jm,n and K^ n'- the current Im,n 
is obtained by substituting the expression for Qm,n ( 12. 39^ into (12.30p and integrating 
with respect to time. Owing to the complexity of the expression for r given by (I2.38p . 
the left-hand side of (I2.30p can not be integrated with respect to time. However, 
the variable r varies more slowly in time than $. Hence, integration by parts (using 
U\t)g{et)<lt = [f{t)g{et)]-Ejf{t)g'{et)dt) gives, to leading-order, 

2ea 2ea 
Im,n ~ [1 - cos(2/i;)] sin $ sech(/3r) sai{2k) cos $ sech(/3r), (2.40) 

UJ UJ 

where we have taken the constant of integration to be zero, and 1] ~ to leading order. 
Similarly, substituting for Qm,n in equations (12.3ip and ( 12. 32^ and integrating, we find 

2ea 2ea 
Jm,n ~ [1 - cos(A; - //i)] sin$sech(/3r) sin(A; - //i) cos $ sech(/3r), (2.41) 

UJ UJ 

2s(x 2scy 
Km,n ~ [1 - cos(A; + Ih)] sin $ sech(/?r) sin(fc + Ih) cos $ sech(/3r). (2.42) 

UJ UJ 

Substituting these expressions into ( 12.291) . we obtain 

2e^a^ 

Eq r^Yl — ^os^ ^ sech^(/3r) 

m,n ^0 

+ —sech^ipr) {[(1 - cos(2A;)) sin $ — sin(2fc) cos $] 

UJ ^ 

+ [(1 - cos(A; - Ih)) sin $ - sin(A; - Ih) cos 

+ [(1 — cos{k + Ih)) sin ^ -sm{k + lh)cos^f] . (2.43) 

We replace the sum by an integral since the variables Z = e{m — ut) and W = e{hn — vt) 
vary slowly with m and n. To simplify the resulting integral, we approximate the terms 
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cos^$, sin^ $ and sin $ cos $ by their average values of |, | and respectively. Hence 
fl2.43p becomes 

^o~E^^sech^W- (2.44) 

m,n ^0 

From the definition of r (12.381) we note the function sech(/?r) is not in general 
radially symmetric in m, n. Hence we work in ?7)-space to facilitate evaluation of the 
double integral which approximates the double sum in f l2.44p . Evaluating the Jacobian 
associated with the transformation from (m, n) to 77) coordinates, we find 

Eo ~ -^^^AD,D2-Dl II sechi/3^e + v') d^dr/. (2.45) 

Evaluating this and substituting for a and /? in terms oi D = 3/2Lj{k,l) and B = 
3buj{k,l)/2 (see the appendix) we find 

47rlog2(21og2 + 1) 1 

It is evident from (12.461) that the leading-order energy Eq is independent of the breather 
amplitude, again confirming the existence of a minimum energy of moving breathers in 
the two-dimensional HETL with symmetric potential. However, the threshold energy 
does depend upon the wavenumbers k and / hence moving breathers have a different 
threshold energy. A plot of the expression (]2.46l) is shown in Figure [7l We see 
that Eq, given by (12.460 . is strictly positive in the region of ellipticity V, and is 
maximised (attaining the same value) at each of the points corresponding to wavevectors 
{ki,...,k6}, that is, at the points which correspond to stationary breathers. The 
threshold energy (12.461) decays to zero towards the boundary of the elliptic domain V 
(see Figure [5|). Hence, as for breathers in the square lattice, the energy threshold for 
moving breathers is lower than that for stationary breathers. The threshold becomes 
arbitrarily small at the boundary of the domain of ellipticity. 

2. 7. Lattices with an asymmetric potential 

In this section, we consider the more general scenario where the potential T '{Q) 
is asymmetric. In this case, the terms a and c in (12. 9p and (12.101) are not both zero. 
If ( 12.2ip is to be reduced to an NLS equation in F, both Go and G2 must be found in 
terms of F. Whilst G2 is given by the simple algebraic equation (I2.20p : in order to find 
Go, the partial differential equation ( 12.23p must be solved. We assume that Go travels at 
the same velocity as F, that is, Go{X,Y,t,T) = Gq{Z,W,T). Eliminating Gom (12.231) 
becomes 

{u^ -Q)Gozz + {v^ -2h^)Goww + 2uvGozw = Qa\F\lz + 2ae\F\lywi2-'^7) 

For general k and /, it is difficult to solve for Go explicitly. However, for any of the 
wavevectors {ki,...,k6}, the velocities u and v become zero, and so (12.471) becomes 
V^Go = — aV^|Fp, where the operator V^, defined by 'V^z,w) = ^1 + ^w^ equivalent 
to Vfx,Y) ^ '^x ~^ '^Y- It follows that Go = — a|Fp for each of these wavevectors. 
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Equations (12.18p -( l2.23h are the same whichever of the wavevectors {ki, . . . , ke} is 
used. Thus, fl2.18p gives = 3, and (12.201) implies G2 = aF'^/3. Substituting these 
expressions for Gq and G2 in (I2.2ip gives the following NLS equation for asymmetric 
potentials for any of the wavevectors {ki, . . . , ke} 

2icuFr + aV^F + uj\3b - fa^)\F\'^F = 0, (2.48) 

The anomalous dispersion regime thus corresponds to 6 > lOa^/9. Soliton solutions are 
derived in the appendix; (I2.48p corresponds to = 1/2 and B = (96 — 10a^)/2, and 
also = + y^. Substituting the solution for F into the lattice ansatz fl2.17p . along 
with the known expressions for Go and G2 gives the following second-order formula for 
stationary breathers in lattices with an asymmetric potential 

Qm,nif) = leacoslkm + Ihn + {oj + e'^\)t] sech(/3r) 

+lae^ a^sech^{(3r)[cos[2km + 2lhn + {2uj + 2e^\)t\ - 3} + O(e^), (2.49) 

where a and [3 are determined in the appendix, r = VA^+T^, and X = em and 
Y = ehn in terms of the original discrete variables m and n. 

2.7.1. Breather energy. We calculate the leading-order energy of stationary 
breathers in lattices with an asymmetric potential. In this case, stationary breathers 
are given by (12. 49^ . from which we take only the leading order term 

Qm,n{t) ~ 2ea cos sech(/5r) , (2.50) 

where 6 = km + Ihn + Vtt is the phase of the carrier wave with (/c, /) corresponding to 
one of {ki, . . . , kg}, Vt = uo + is the breather frequency, and w = ?>. 
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The currents /m,n, Jm,n and Km,n are obtained by substituting (12.501) into equations 
fl2.30p - fl2.32p and integrating with respect to time, taking the constant of integration to 
be zero. Thus, to leading order, we find 

lm,n ~ sinBsech(pr) cos B sech(pr), (2.51) 

UJ uo 

3ea 1 / /I \ \/^sa ^ f n \ 

Jm,n ~ sm B sech(pr) H cos B sech(pr), (2.52) 

OJ UJ 

Sea , / n \ V^^ce ^ 1 / /I \ 
Km,n ~ sm B sech(pr) H cos B sech(pr). (2.53) 

CO UJ 

Substituting these expressions into f l2.29p gives the leading-order energy 
^0 ~ E cos' B sech2(/3r) 



m,n 

,2 



Lh 



3h^ sin' e + 2h sin B cos B + 3 cos' B sech'(/3r). (2.54) 



We approximate the term in square brackets by taking the average values of cos'B, 
sin'B, sin B cos B as |, | and respectively. 

We replace the double sum by an integral giving 

- // sech'(/3VX2 + y2) dXdy. (2.55) 









HQ 

This can be evaluated to 

47rlog2a'_ Stt log 2(2 log 2 + 1) 
° /Z\r^~ Co/i(96-10a2)(41og2-l)' 
Again, this estimate for energy is independent of the breather amplitude A, 
demonstrating the energy threshold properties of the two-dimensional HETL; namely, 
the activation energy required to create a breather in the HETL is an 0{1) quantity, 
irrespective of its amplitude {2ea ^ 1). 



3. Higher-order asymptotic analysis 

The cubic NLS equation exhibits blow-up, which cannot occur in a discrete system 
since energy is conserved, and even if all the energy were localised at a single site, the 
amplitude would still be finite. The two-dimensional cubic NLS equation description of 
the breather envelope is lacking. Higher-order dispersive and nonlinear effects play 
an important role in the dynamics of such discrete systems and therefore must be 
incorporated. In this section we extend our analysis of the lattice equations (12.10p to 
fifth-order, and derive a generalised NLS equation which includes higher-order dispersive 
and nonlinear terms. It then remains to determine whether this generalised NLS 
equation supports stable soliton solutions. 

Due to the complexity of a fifth-order analysis of a general asymmetric potential, 
we consider only lattices with a symmetric potential, that is, those for which a = c = 
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in (12. 101) . Since no second or fourth harmonic terms are generated by the nonlinearity, 
we use a much simpler ansatz, namely 

Qn^,n(^) = e^'^F{X, Y, r, T) + e''e^'^H^{X, y, r, T) + ■ ■ ■ + c.c, (3.1) 

where the phase ifj = km + Ihn + ut. In this case, in addition to the equations fl2.24p - 
(I226D, (12341 ) and ^M>, we also have 

0(5^^): 

Ftt = I [8 cos(2fc) + cos k cos{lh)] Fxxxx + cos k cos{lh)FxxYY 

+ I /i^ cos /c cos (//i)FYyyy — ^h siu k sm{lh)FxxxY — sin A; sin(//;,)Fxyyy 
+ 66 [2 cos(2A;) + cos k cos{lh)] {\F\'^F)xx + Qbh^ cos k cos{lh){\F\^F)YY 

— 126 sin ksm{lh) hhFFxFy + 2hFFYFx + 2hFFxFY + hF^FxY + ^HFFFxy] 

- Sbuj^F^H^ - 10uj^d\F\^F. (3.2) 

We simplify this by restricting attention to stationary breathers. Accordingly, only one 
extra timescale, T = eH, is required, and we fix the wavenumbers k and / to correspond 
to one of the wavevectors {ki, . . . , kg}. Hence we obtain the equations: 



2icjFT + 3V^F + 36cj2|FpF = 0, (3.3) 



Ftt = -fV^F - 9bVWF\^F) - 10uj^d\F\^F, (3.4) 

in addition to = 3 (from the 0{ee^'^) equation) and = (from C(5^e^'^)). A 
consequence of the hexagonal symmetry of the HETL is that all differentials on the 
right-hand side of (13. 4p are isotropic. 

In order to obtain a generalised NLS equation, we combine the higher-order equation 
(13. 4p with the cubic two-dimensional NLS equation (]3.3I) . The Ftt term on the left-hand 
side of ( 13. 4p is eliminated by differentiating ( 13. 3p with respect to T and substituting for 
Ftt into (13.41) . The resulting expression is 

6iFT + 3V^F + 27b\F\^F + —V^F + — (40d - 21}?)\F\^F 

+1!|^V2(|F|2f) - :^\F\^V^F - ^-^F^V^F = 0, (3.5) 

which is isotropic. 

To the best of our knowledge, this perturbed form of the NLS equation has not 
been studied in the literature before; although it is similar to and slightly simpler than 
the corresponding equation derived for the square lattice It is also similar to the 
perturbed NLS equation considered by Davydova et al namely 

iFt + DV'^F + B\F\^F + PV^F + K\F\'^F = 0, (3.6) 
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which is known to have stable soliton solutions. The anomalous dispersion case is 
BD > 0, which in our equation (|3.5I) corresponds to 6 > 0; Davydova's criteria for 
soliton existence is PK > which implies AOd > 276^. Hence, it is in this parameter 
regime that we seek breathers solutions of the HETL. The numerical results presented in 
Section [4] show that long-lived breather solutions are supported by the two-dimensional 
hexagonal FPU lattice, suggesting that the additional perturbing terms in (I5L5I) do not 
destabilise the Townes soliton. In fact, numerical simulations (presented in [l| but not 
here) of the case AOd < 276^ show the breather mode to be long-lived, suggesting that 
the terms on the second line of (13. Sp are stabilising. 

We have been unable to find a variational formulation of (13. Sp . and are therefore 



unable to use the methods of Davydova et al. [J| and Kuznetsov et al. [13|. Alternative 
possible methods include the modulation theory of Fibich and Papanicolaou 0, 0|- 

4. Numerical results 

4.1. Preliminaries 

Using a fourth-order Runge-Kutta scheme, we solve the equations 
dQ 

m.n r-> 

dt 
d /? 

= [5] + 5] + 51) [g„,„ + aQl^^ + hQl^^ + cQ^,„ + dQl^^^ , (4.1) 

numerically. Introducing the variable Rm,n converts the system of second-order ordinary 
differential equations (I2.10p to an equivalent system of first-order differential equations. 

We present the results of simulations for a range of parameter values. From (I2.26p . 
the velocity of the envelope {u,v) depends upon the wavevector k = [/c,/]"^, and hence 
we obtain moving breathers by choosing {k, I) G V. In Figure [HI we show the points 
in Pi C P for which we solve the lattice equations (14. ip numerically. These points 
correspond to the wavevectors ki = [7r/3, 7r//i]-^, ka = [1.4, 7r//i]^, kb = [0.79,1.7324]^ 
and kc = [0.8,1.9987]^. Selecting (/c, /) too near to the boundary of V results in a 
sharply elongated breather, which are difficult to simulate as they require a large domain. 
Wavevectors close to any of {ki, . . . , kg}, lead to breather modes with small speeds, a 
check of breather velocity would then require a long-time simulation. In practice we 
have chosen wavenumbers which do not lead to a severely elongated breather envelope, 
and yet have velocities which result in observable displacements over reasonable times. 
We present simulations of stationary breathers in systems with asymmetric potentials 
(a, c 7^ 0) in Section 14.81 and, in Sections 14. 4114. 71 simulations of stationary and moving 
breathers in lattices with symmetric potentials (that is, a = c = in (14. ip ). 

4-2. Initial data and boundary conditions 

We generate initial data by using the analytic expressions for breather solutions 
derived in Section 12.31 The formulae for Qm,n and Rm,n are found in terms of the 
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Figure 8: Wavevectors in Pi C P for which breathers are simulated. 

original discrete variables m and n, and then shifted horizontally and vertically so that 
initially the breather lies at the centre of the lattice. We impose periodic boundary 
conditions for the lattice in both horizontal and vertical directions, converting the 
two-dimensional arrangement illustrated in Figure [U into a two-torus. In long-time 
simulations, moving breathers which approach an edge of the lattice reappear from 
the opposite edge. We select the site (1, 1) to lie at the bottom left-hand corner of 
the arrangement as illustrated in Figure [9l Sites along the boundaries and corners 
are missing between one and four neighbours. We introduce fictitious sites along the 
boundaries and corners where necessary to effect periodic boundary conditions. 

An illustration of a small finite lattice is shown in Figure [9l The dots represent 
capacitors located at lattice sites, and the lines represent some of the inter-connecting 
inductors. From equations (j4.1l) . the charge Qm,n stored on each capacitor depends 
upon the charge stored on the capacitors located at its six neighbouring sites, two in 
each of the directions Ci, ej and Sk. For the sake of clarity, the inductors connecting 
the capacitors at the center of each hexagon to its six nearest neighbours are not shown 
(see Figure [I]). It is not necessary that the lattice should be "square," meaning that 
the lattice could comprise M x N lattice sites, with M ^ N. However, typically, we 
consider M = as in Figure [9l since the numerical routines are simpler to encode when 
the lattice is square, and in the examples below, we consider lattices with < 50. 

4-3. Numerical computation of breather energy 

Since the HETL is lossless the total energy is conserved and can be used as a check 
of the accuracy of our numerical scheme. We compute the leading-order energy by 
expressing the summand e^^„ f l2.29p in terms of the output variables of the numerical 
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Figure 9: Periodic boundary conditions for the two-dimensional HETL. 



routine, namely Qm.^n and Rm,n- The first term of is simply dependent on Qm,n- 
It remains to find the currents Im,n, Jm,n and Km,n in terms of Qm,n and Rm,n- The 
details of this calculation depend upon whether the interaction potential is symmetric 
or asymmetric. As a check of the numerical scheme we verify that the sum is conserved 
and that it agrees with the asymptotic estimate (12.46p . 

4.3.1. Lattices with a symmetric potential. First, we find the current Im,n in terms of 
Qm,n and Rm,n. Differentiating the leading-order expression for the breather given by 
fl2.39p . we have (retaining leading-order terms only) 

Qm,n = Rm,n ~ —2eaujsm^ sech(/3r), (4.2) 

where a, /3 and are as defined in Section 12.61 Comparing the analytic expression 
(I2.40p for the current Im,n with the expressions for Qm,n (12.39^ and Rm,n (I4.2p . we find 

_ cos(2A:) - 1 sin(2fc) 

Similarly, comparing the expressions f l2.4ip and f l2.42p for the currents Jm,n and Km,n 
respectively with equations (12.391) and (14.21) . we find 

_ cos(fc -lh)-l sm{k-lh) 

'Jm,n 9 J^m,n '^m,ny \ ) 

_ cosjk + lh)-l sm{k + lh) 

Substituting these expressions into (12. 29^ . gives an expression for Eq which does not 
simplify, so we do not reproduce it here. 
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4.3.2. Lattices with an asymmetric potential. Differentiating the leading-order 
expression (I2.50p for the charge Qm.n in lattices with an asymmetric potential gives 

Qm,n = Rm,n ~ —2eaujsmQ sech(/3r), (4.6) 

where 6 = 7rm/3 + vrn + f2t with = 3. Using fl2.50p and (14.61) . the current Im,n given 
by (I2.5ip can be expressed in terms of Qm,n and Rm,n as 



2uj 



UJ 



~l~ Qm,r, 



Similarly, from fl2.52p - fl2.53p . we have 



J -K 

'Jm,n J^m.n ~ „ 

2uj 



m,n 



UJ 



Qm,n 



(4.7) 



(4.8) 



These are equivalent to substituting k = ir/S and I = n/h into f l4.3p -f l475l ). Inserting 
the expressions for Im,n, Jm,n and Km,n ( fl4.7p - fl4.8p ) into fl2.29p yields 

1 



eo = y: 



m,n 



72Cn 



^^Qm n ~^ '^R-m n '^^/^Q m,nRm,r. 



(4.9) 



4.3.3. Effective breather width. Numerical computation of the energy as described 
above allows us to check that the total lattice energy is conserved, but does not indicate 
whether a breather changes shape over time. To remedy this, we define breather widths 
in the m and n directions, the sum of which we denote Wbr, where 



br 



(4.10) 



Eq Eq \EqJ 
and rio, roi, and ro2 are defined by 

The variation in Wbr over time gives a measure of the distortion suffered by a breather. 



2 



'l"20 — Y.r. 
^02 X/m.n ^ ^ (^m,n. 



4.4. Stationary breather in a lattice with symmetric potential 

We investigate stationary breather solutions in the anomalous dispersion regime 
which corresponds to 6 > (see f|3.5p ). We set k = n/S and / = n/h, corresponding to 
ki in Figure m so that the velocities u and v are zero. Davydova's result implies that we 
expect to find stable soliton solutions when PK > 0. From (13.60 . P = s'^/2, hence this 
inequality implies iOd > 276^. Hence we choose b = d = 1 as well as = 30, e = 0.2 
and A = 1. Using the technique outlined in the appendix we calculate a = 1.0212, 
P = 1.8670. The breather frequency is + e'^X = 3.040, and therefore the period of 
oscillation is T = 2.0668. 

As is the case for all our simulations, the initial profile of the breather is located at 
the centre of the lattice, as illustrated in Figure [10] (a). At t = 0, the breather energy is 
Eq = 0.7606, and Wbr = 3.74. The asymptotic estimate for Eq given by f l2.46p is 0.7523, 
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which is only 1% different from the numerically obtained value. In Figure fTOl(c) we show 
the breather after 30 oscillations. Plots of the local energy em,n at t = and t = SOT are 
presented in Figures fTOlfb) and [10 (d). After thirty oscillations the breather has shed a 
small amount of energy, which is manifested as small amplitude radiation throughout 
the lattice. Accordingly, the breather appears a little distorted in shape compared to its 
initial profile. In particular, att = SOT, we find Wbr = 4.21. The energy of the breather 
at t = SOT is 0.7087, giving AEq/Eq = -0.0682, 




(c) Profile at t = SOT ^ 62.82. (d) Plot of e„,„, Ea = 0.7087. 

Figure 10: Stationary breather in a lattice with symmetric potential, see Section [44] 
for details. 

4-5. Breather moving along a lattice direction ("^ = 0°) 

In Figure [TT], we show a simulation of a breather moving along a lattice direction 
parallel to the m-axis, that is, \E' = 0°. We have chosen k = ka = [1.4, tt/Zi]-^ so that 
u = 0.4445 and v = and thus ^ = tan~^(i;/n) = 0°. The breather frequency is 
uj = 2.9265, and hence the period T = 2.1470. Following the calculation outlined in 
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the appendix, the amplitude and width parameters are a = 1.0339 and (3 = 1.8440; the 
remaining parameters being b = 1, d = 1, N = 30, e = 0.1 and A = 1. 

The initial profile of the breather is shown in Figure fTT](a). and at this time, the 
calculated energy is £"0 = 0.5537, whilst the asymptotic estimate (12.461) is Eq = 0.5522. 
It may be observed that the breather is not radially symmetric. In fact, it is slightly 
elongated in the direction parallel to the m-axis, that is, parallel to the direction of 
motion. This is because the point corresponding to the wavevector ka is near the 
boundary of the region of ellipticity Vi as illustrated in Figure O 

Figure[Tl](b) shows the breather at the later time of t = 21.58T = 46.3513, at which 
time we find Eq = 0.5562. By this time, the breather has reached the right-hand edge 
of the lattice and, owing to periodic boundary conditions, it reemerges from the left- 
hand side as shown in Figure fTTlfc). The breather remains localised, without spreading, 
though it leaves behind a small amount of energy in its path (visible in Figures fTT](b) 
and [TTlfc)). The computed energy changes much less than for the previous simulation 
(here, AEq/Eq = 0.0287). 

The velocity of the breather can be measured from a plot of the energy e^.n, which 
is shown in Figure [ll](d), viewed from directly above the plane of the lattice. From this 
plot, we note that the breather has travelled 42 units at an average speed of 0.42 units 
per second, which is 5.5% lower than our predicted speed of 0.4445 units per second. 
This is consistent with a small amount of energy being shed as the system transforms 
from our approximated initial conditions into the precise shape of the breather. 

4.6. Breather moving at = 210° 

We now show that it is possible to simulate breathers moving in directions other 
than a lattice direction (that is, \Ef 7^ 0°). We set k = 0.79 and / = 1.7324, which 
corresponds to kb in Figure El It may be verified that u = —0.1999 and v = —0.1154 
units per second, leading to = 210° as required. Although not a lattice vector, 
this direction is an axis of symmetry of the lattice. For the wavevector kb, we have 
u = 2.9675 and hence T = 2.1174, the remaining parameters being b = 1, d = 1, 
N = 30, e = 0.1 and A = 1. The variational parameters a and (3 are 1.0267 and 1.8568 
respectively. The breather is shown at times t = 25, 50, 75 and 100 seconds in Figure [121 
Clearly, the breather does not deform significantly as it travels, nor does it radiate much 
energy. The initial energy is computed to be Eq = 0.6306, and at t = 100, the energy 
is 0.6275, a loss of 0.5%. The asymptotic estimate for the energy is Eq = 0.6184 — 2% 
different from the numerically computed value. The motion of the breather is charted 
in Table [H the final measured values for the velocities u and v give an average speed 
of 0.2152 units per second, 6.8% below the predicted speed of 0.2308 units per second. 
The angle of travel is almost identical to the expected value of \& = 210°. 
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(c) Profile at t = 46.59T, Eq = 0.5696. (d) Plot of e™,„. 

Figure 11: Breather moving along a lattice direction, ^' = 0°, see Section [4?5l for 
details. 



Time 

(s) 


Horizontal 
displacement 


Vertical 
displacement 


Average horizontal 
velocity (units s"""^) 


Average vertical 
velocity (units s"""^) 


tan * 




25 


-3.5 


-2 


-0.14 


-0.08 


0.5714 


209.74° 


50 


-8 


-4.5 


-0.16 


-0.09 


0.5625 


209.36° 


75 


-13 


-7.5 


-0.1733 


-0.1 


0.5770 


209.99° 


100 


-18.5 


-11 


-0.185 


-0.11 


0.5946 


210.74° 



Table 1: Summary of breather motion = 210°). 



4-7. Breather moving at = 130° 

We have presented simulations of breathers v^hich move along axes of symmetry of 
the lattice. In Section 12.6.11 from the results of our asymptotic analysis, we found that 
breather solutions could be constructed for any direction of travel. In this section, we test 
the mobility of breathers along directions which do not correspond to axes of symmetry 
of the lattice. We have successfully propagated breathers in a range of directions, and 
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here we show only one such breather, moving at an angle \i/ = 130°. We set k = 0.8 and 
/ = 1.9987, from which we find u = -0.2160, v = 0.2575, * = 130°, u = 2.9502 and 
the period is T = 2.1298. The remaining parameters are b = 1, d = 1, N = 30, e = 0.1, 
A = 1, a = 1.0297 and (3 = 1.8514. From Figure [131 Even at 120 seconds, the breather 
remains localised without suffering appreciable degradation nor is much radiation left 
behind in its wake. Initially the energy Eq is computed as 0.5977, only 0.2% different 
from the asymptotic estimate of Eq = 0.5965. After 120 seconds, Eq is 0.5897, a change 
of 1.3%. To find the velocity of the breather, we have recorded its motion and presented 
relevant data in Table [2l From (12.261) . the average speed is predicted to be 0.3254 units 
per second; whilst our numerical simulation shows an average speed of 0.3361 units per 
second, 3.2% lower than the expected speed. The direction of motion of the breather is 
predicted accurately. 
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" m " m 

(c) Profile at t ^ 90, Eq = 0.5939. (d) Profile at t = 120, Ea = 0.5897. 

Figure 13: Breather moving at ^' = 130°, see Section l47l for details. 



Time 
(s) 


Horizontal 
displacement 


Vertical 
displacement 


Average horizontal 
velocity (units s""*^) 


Average vertical 
velocity (units s~^) 


tan ^ 




30 


-6 


7 


-0.2 


0.2333 


-1.1667 


130.60° 


60 


-11 


13.5 


-0.1833 


0.225 


-1.2273 


129.17° 


90 


-17.5 


21.5 


-0.1944 


0.2389 


-1.2286 


129.14° 


120 


-25 


30 


-0.2083 


0.25 


-1.2 


129.81° 



Table 2: Summary of breather motion ("if = 130°) see Section [471 for details. 



4-8. Stationary breather in a lattice with asymmetric potential 

In Sections 14 .41 - 14 .71 we have shown simulations of lattices with a symmetric 
potential, that is, with a and c not necessarily zero in ( I4.ip . We now consider the 
more general case for which the interaction potential is asymmetric (a 7^ 7^ c). 
To illustrate a stationary breather, we generate initial data using (12.491) . with the 
wavevector k = ki = [tt/S, tt/Zi]-^. As discussed in Section [2771 anomalous dispersion 
corresponds to b > lOa^/9, hence we choose a = 1, b = 2.5, c = and d = 1. The 
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remaining parameter values are = 30, e = 0.1, A = 1, hence = 3, T = 2.0668, 
a = 0.8665 and p = 1.8670. 

The breather is shown in Figure [H] after 10, 30 and 40 full oscillations. Initially, we 
find the breather's width and energy are Wbr = 7.46 and Eq = 0.5429, the latter being 
only 0.26% different from the asymptotic estimate f l2.56p of 0.5416. The accompanying 
plots of Cm.n demonstrate clearly that the breather preserves its form and remains 
localised, even after 80 seconds. At t = 40T, we find that Wbr = 6.83, a narrowing 
of the breather by 8%. Also, at t = 40T, the numerically computed value of the energy 
Eq = 0.5371, showing that the energy does not fluctuate significantly. 



5. Discussion 



In this paper, we have found approximations to discrete breathers in a two- 
dimensional hexagonal FPU lattice with a scalar-valued function at each node. We 
have shown that the lattice equations can be reduced to a cubic NLS equation for 
two special cases. In Section 12.61 we obtained moving breathers when the interaction 
potential is symmetric; in this case an ellipticity criterion for the wavevector is found. 
In the appendix we summarise our technique for approximating soliton solutions of the 
two-dimensional NLS equation. 

In Section 12.71 we considered lattices with an asymmetric potential, in which case 
a reduction to NLS could only be performed for stationary breathers. The theoretical 
methods employed here are similar to those used on the square lattice of however, 
several important differences emerge in the course of the analysis. We found that 
the anomalous dispersion regime corresponds to b > 4a^/3 in the square lattice and 
b > lOa^/9 in the hexagonal. Furthermore we find that stationary breathers involve the 
generation of second harmonics in the hexagonal lattice {G2 = aF^/3), whereas they are 
suppressed in the square lattice {G2 = 0). In both lattices, the quadratic nonlinearity 
generates a small amplitude slowly-varying mode (Go = 

For symmetric interactions we found an associated ellipticity constraint: moving 
breathers only occur for certain wavevectors, this means that all breathers have (i) 
a relatively high frequency, and (ii) a maximum speed, which depends on the wave 
vector and hence on the direction of travel, (iii) a threshold energy which also depends 
on wavevector and hence is related to speed and direction of travel. We find that 
the threshold energy is lower for breathers at the edge of the domain of ellipticity 
in wavevector space. These breathers have lower frequencies and faster speeds. The 
fastest-moving breathers are restricted to moving along lattice directions. 

We also presented asymptotic estimates for the breather energy in Sections 12.61 
and 12. 7[ As expected, we found a minimum energy required to create breathers in 
the hexagonal lattice. The threshold energy for moving breathers is smaller than that 
required for stationary breathers and becomes vanishingly small at the boundary of the 
domain of ellipticity. 

In Section m we extended the small amplitude expansion to fifth-order and derived a 
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Figure 14: Stationary breather in a lattice with an asymmetric potential, see Section 
[48] for details. 
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higher-order equation which more correctly describes the shape and stability properties 
of the breather envelope. We obtained a generalised NLS equation (13.51) with a variety of 
perturbation terms, some of which are known to be stabilising. This equation is slightly 
simpler than the corresponding equation obtained for the square lattice, namely (3.11) 
of In particular, the higher-order dispersive terms in (13.51) are isotropic, reflecting 
the hexagonal rotational symmetry of the lattice. For stationary breathers in the case 
of a symmetric potential we flnd that the cubic nonlinearity does not give rise to third 
harmonics, that is, ifa = in the hexagonal lattice, in contrast to = bF^/8 in the 
square lattice. 

In Section m we illustrated these breather modes, showing that both stationary and 
moving breathers are long-lived. The breather proflles change little over time, a small 
amount of energy is shed due to the initial conditions being only approximate. We have 
successfully propagated long-lived breathers moving in directions which are not axes of 
symmetry of the lattice (for instance, \Ef = 130°), suggesting that there is no absolute 
restriction upon the direction of travel; this is in contrast to the observations reported 
by Marin et al 15, [3| for mechanical (two-component) two-dimensional lattices, who 



could only flnd breathers which travelled along axes of symmetry of the lattice. 
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Appendix A. Approximation to Townes soliton 

Since analytic formulae for Townes solitons are unavailable, we use the Rayleigh- 
Ritz method to flnd time-harmonic radially symmetric solutions of 

iFT + DV^F + B\F\'^F = 0, (Al) 

of the form F(x, T) = e''^"^0(r) where r = |x| = + v"^- The function satisfles 

-A0+DVV + fi0^ = O, (A2) 

where = 9| + 9^. Equation (1A2[) arises from a variational derivative of 

= J j\m' + i^lV^r - \B\ct>\'d?v. (A3) 

Using a trial solution of the form = Q!sech(/3r), where a and {3 are undetermined 
parameters, we flnd 

L'(l + 21n2) 2 5(41n2-l)a4 Aln2 ^ 

a and (3 are determined by seeking stationary points of the action £, namely d£/da = 
d£/dp = 0. Hence we flnd 

I 12Aln2 I 6Aln2 

"-V5(41n2-1)' V/^(21n2 + l)' ^^^^ 
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and so our approximation to the Townes soliton solution of ( lAll ) is 

F = J ^""^ , exp(iAT) sech ( ■ (A6) 

References 

[1] lA Butt. PhD Thesis. University of Nottingham, to be published at 

|http : / / etheses . nott ingham . ac . uk / (2006). 
[2] lA Butt and JAD Wattis. Discrete breathers in a two-dimensional Fermi-Pasta-Ulam lattice. J. 

Phys. A: Math. Gen., 39, 4955, (2006). 
[3] lA Butt and JAD Wattis. Asymptotic analysis of combined breather-kink modes in a Fermi-Pasta- 
Ulam chain, submitted, (2006). 
[4] TA Davydova, AI Yakimenko, and YA Zaliznyak. Two-dimensional solitons and vortices in normal 

and anomalous dispersive media. Phys. Rev. E, 67, 026402, (2003). 
[5] E Fermi, J Pasta, and S Ulam. Studies of nonlinear problems. Los Alamos Scientific Report, 1940, 

originally unpublished. Later published in Lectures in Applied Mathematics, 15, 143, (1974). 
[6] G Fibich and G Papanicolaou. A modulation method for self-focusing in the perturbed critical 

nonlinear Schrodinger equation. Phys. Lett. A, 239, 167, (1998). 
[7] G Fibich and G Papanicolaou. Self-focusing in the perturbed and unperturbed nonlinear 

Schrodinger equation in critical dimension. SIAM J. Appl. Math., 60, 183, (1999). 
[8] S Flach, K Kladko, and RS MacKay. Energy thresholds for discrete breathers in one-, two-, and 

three-dimensional lattices. Phys. Rev. Lett., 78, 1207, (1997). 
[9] S Flach, K Kladko, and CR Willis. Localized excitations in two-dimensional lattices. Phys. Rev. 

E, 50, 2293, (1994). 

[10] J GiannouHs and A Mielke. The nonhnear Schrodinger equation as a macroscopic limit for an 

oscillator chain with cubic nonlinearities. Nonlinearity, 17, 551, (2004). 
[11] J Giannoulis and A Mielke. Dispersive evolution of pulses in oscillator chains with general 

interaction potentials. Disc. Cont. Dyn. Sys. B, 6, 493, (2006). 
[12] M Kastner. Energy thresholds for discrete breathers. Phys. Rev. Lett., 92, 104301, (2004). 
[13] EA Kuznetsov, AM Rubenchik and VE Zakharov. Soliton stability in plasmas and hydrodynamics. 

Phys Rep, 142, 103, (1986). 
[14] RS MacKay and S Aubry. Proof of existence of breathers for time-reversible or Hamiltonian 

networks of weakly coupled oscillators. Nonlinearity, 7, 1623, (1984). 
[15] JL Marin, JC Eilbeck, and FM Russell. LocaHsed moving breathers in a 2D hexagonal lattice. 

Phys. Lett. A, 248, 225, (1998). 
[16] JL Marin, JC Eilbeck, and FM Russell. 2D breathers and applications (Nonlinear Science at the 

Dawn of the 21st Century). Eds: PL Christiansen and MP Soerensen. Springer, Berlin, (2000). 
[17] JL Marin, FM Russell, , and JC Eilbeck. Breathers in cuprate-like lattices. Phys. Lett. A, 281, 

21, (2001). 

[18] M Remoissenet. Low-amplitude breather and envelope solitons in quasi-one-dimensional physical 

models. Phys. Rev. B, 33, 2386, (1985). 
[19] MI Weinstein. Excitation thresholds for nonlinear localized modes on lattices. Nonlinearity, 12, 

673, (1999). 



